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We study the viscosity and thermodynamic properties of QGP at RHIC by em- 
. ploying the recently extracted equilibrium distribution functions from two hot QCD 

equations of state of 0(g 5 ) and 0(g 6 ln(l/g)) respectively. After obtaining the tem- 
^ | ■ perature dependence of energy density, and entropy density, we focus our attention 

on the determination of shear viscosity for a rapidly expanding interacting plasma, 
as a function of temperature. We find that interactions significantly decrease the 
shear viscosity. They decrease the viscosity to entropy density ratio, ij/5 as well. 
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Recent experimental results from RHICpH] reveal that the QGP produced in heavy ion 

collisions behaves like an almost perfect fluid with very low viscosity Q, fl Q. In belying 
c3 ■ 

the earlier expectations that the deconfined phase would show nearly ideal behavior at 
temperatures close to T c , the results from the flow measurements signal that the deconfined 



phase is strongly interacting. Lattice studies 



5| also predict that the equation of state for 



QGP is about 10 percent away from ideal EOS even at T ~ 4T C . Therefore, studies based 
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on the ideal EOS in the computations are certainly inadequate to address the physics in this 
domain of QCD. 

In an attempt to understand the flow measurement results, a variety of techniques have 
been employed, including the employment of AdS/CFT correspondence. A strong result 
from these studies is a lower bound on the viscosity to entropy density ratio given by rj/S < 
j- ~ 0.08 {2]. Yet another approach is to study a strongly coupled classical plasma jy] which 
again yields a value close to the lower bound mentioned above. There are two lattice results, 
one by Meyer [3] who obtains a value 0.13 for rj/S in pure gauge theory at T = 1.56T C . The 



other result, also in pure aguge theory is due to Nakamura 



8| who finds that rj/S < 1. 



Recall that the analyses 

@,T3 based on v 2 measurements |H] arrive at values that vary 
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ies 0, Q, Q, 



from 0.08 — 0.2. Results of other studies [3|, |4|, [12J, [13] also yield numbers in the same range. 
Interestingly, Asakawa et al H] find that the ratio can take a value which is smaller 

than the lower bound set by AdS/CFT studies, depending on the value assigned to the 
transport parameter q^. 

It is noteworthy that while the perfect liquid picture implies a strongly interacting QGP 
(sQGP), most of the estimates made above employ a perturbation about the ideal distri- 
bution for quarks and gluons. While this could simplify matters, it is important to find 
out what an EOS which includes the interactions would predict for rj/S. In doing so, one 
could also gain an appreciation of the nature of interactions in sQGP. In th is pape r, we 



study the predictions of EOS which are based on improved perturbative QCD 17|, ll8(. The 
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ications of lattice EOS will be taken up in a subsequent paper. Our study is based on 



201 ] where it has been shown how the EOS may be adapted to study the properties of 



QGP in heavy ion collisions. It is worth mentioning that the above mentioned works have 
demonstrated the viability of the EOS since they yield reasonable values for the dissociation 
temperatures for J/\l/ and T. 

The paper is organized as follows. In Section (II), we introduce the EOS based on pQCD 

nn 

and review the work contained in pj3|, 12DJ • We proceed to determine the the temperature de- 
pendence of thermodynamic observables (energy density and entropy density) in Section(III). 
In Section (IV), We obtain the expressions for anomalous and collisional contributions to 
the parton shear viscosities, in the presence of interactions. We further study the behavior 
of viscosity to entropy density {i]/S) as a function of temperature for pure gauge theory 
plasma, and compare the results with the ones obtained from an ideal gas distribution. In 



Section(VI), we conclude this work. 



II. HOT QCD EQUATIONS OF STATE AND THEIR QUASIPARTICLE 

DESCRIPTION 

Recently Chandra et. a/fiol. I20I] have considered two EOS based on pQCD, and developed 
a self-consistent method to recast them in terms of non-interacting/weakly interacting quasi 
particles with effective fugacities. Since the method is employed here, we briefly review the 
work. 

The EOS which we label EOS2 18j] is given by 
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The other EOS |l7j, which we call EOS1, is of 0(g 5 ), and is obtained from this equation by 
dropping the last term which has contributions of 0(g 6 (ln(l/g)+5)). The phenomenological 



parameter 5 is introduced in 18| to incorporate the undetermined contributions of 0(g 6 ). 

The above equation of state has several ambiguities, associated with the renormalization 
scale (hms)i the scale parameter At/A^ s which occurs in the expression for the running 
copulping constant a s , and the value of the phenomenological parameter 5. The ambiguity 
associated with (fiMs) has been discussed well in literature and a popular way out is the 



BLM criterion due to Brodsky, Lepage and Mackenzie 



In this criterion, which is chosen 



to make the highest power of Nf vanish in the highest perturbative order, the value of {^ms) 



is allowed to vary between nT and 4nT 



24l | . In this paper we choose the renormalization 



scale /ijv/5 = 2.157rT ph 6.752T25J close to the central value 2nT. One feature of this 
particular choice is that all the contributions due to the logarithms containing fi^ s are very 
small. For the scale parameter At, we follow Huang and Lissia 22J and set Kt/^ms = 
exp(7£ + l/22)/47r rs 0.148, since with this choice, they find among other things that the 
coupling g 2 (T) is optimal for lattice perturbative calculations. The same value has also been 
emfployed by others, see e.g. [18j, [25|. Finally, we set Aj^s = T c , which is close to the value 
0.87T C found by Gupta Q- 

We turn our attention to the phenomenological parameter 5. The optimal value of 5 
depends on the choice of the renormalization scale and the order in which the running 



coupling constant is determined. Blaizot, Iancu and Rebhan 



261 ] find that the optimal value 



is given by 5 = 1/3 if one employs the two loop running coup 



loop running coupling constant yields 5 in the range 0.7-0.9 18j, |20j]. In this paper, we find 



ing, constant while, the one 



that 5 in the range 0.8 to 1.2 yields the best fit with the lattice results. The important point 
here is that once the phenomenological parameter S is fixed by comparing EOS2 with lattice 
EOS, it can be employed to study the properties of QGP. As regards EOS1, we note that 
the matching with the lattice results has been found to be merely qualitative 
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20|. In 



Let us briefly review the underlying idea and the findings of recent two papers 
Ref . [19I] . it has been shown that the interaction effects in EOS1 and EOS2 can be captured 
in terms of the effective fugacities {z g ,z q ) for the quasi gluons and quarks. The effective 
fugacities are determined self-consistently order by order. The mapping has been found to 
be accurate up to about 5% error. Therefore, we expect an error of the same order for all 
the quantities which can directly be derived from the pressure for eg., the energy-density 
and the entropy density. 

Interestingly, the temperature dependence of the screening length which is subsequently 
determined is seen to qualitatively agree with the lattice results of Zantow28(. In Ref. 20], 



the quasi-particle description deve 
of the response function of QGP 



oped in Ref. 19J has been combined with the formulation 



29], and the dissociation temperatures for J/\& and T 



have been estimated. These numbers are again reasonably close to the predictions of other 



theoretical works 



311 ] . This motivates us to further utilize EOS1 and EOS2 to study the 



behavior of thermodynamic quantities such as energy density, entropy density, and most 



importantly, the transport parameters, shear viscosity rj and viscosity to entropy density 
ratio rj/S, for the rapidly expanding plasma. In addressing this, we generalize the recent 
work of Asakawa, Miiller and Bass 14J on the transport properties of interacting QGP. 



III. THERMODYNAMIC OBSERVABLES 



Let us now turn our attention to study the behavior of thermodynamic observables. We 
consider energy density (e) first which brings out the physics of the quasi-particle description 
manifestly. We then determine the entropy density (S), for both EOS1 and EOS2. We 
principally employ the method developed in [ig| ] . 

As mentioned, EOS1 and EOS2 are mapped to the corresponding equilibrium functions 
with the quarks and the gluons possessing effective fugacities: 
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(2) 



where all the interaction effects are captured in the fugacities z g / q = exp/3fi g / q . The form 
of function of temperature is given in [20[. With these distributions, it is straight 

forward to determine the thermodynamic quantities. 



A. The energy-density 

Notwithstanding appearances, the energy of the quasi-gluons and quasi-quarks in not 
merely given by the relation E p = p. Rather, it should be determined from the fun- 
damental thermodynamic relation between the energy density and the partition function, 
e = —dp\n(Z). Substituting for the partion function in terms of quasi-gluons and quasi- 
quarks we obatin, 

6(1/9 = ^I iP + ^ ln ^)^ ( 3 ) 
where v = (u g , u q ) = (2(N 2 — l),4N c Nf). The modified dispersion relation reads, 

E p = p + T 2 d T ln(z g/q ). (4) 

After performing the momentum integration in Eqj3j we obtain the following expression 
for the energy- density: 

Ya = ^6VolyLog[4,z g )}-^6VolyLog[4,-z q )} 
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where 



(5) 



A, = T 2 d T \n(z g )Af g 
A q = T 2 d T ln(z q )N q , 



(6) 



are the contributions from the quasi-gluons and quasi-quarks to the trace anomaly. The 
second term in the dispersion relation Eq.fllJ) may be thus looked upon as the anomalous 
component of the dispersion relation. The quantities Af g and Af q are the quasi-gluon and 
quasi-quark number densities and having the following form, 



v T 3 

M g = -2-r-VolyLog[3,z g 
v T 3 

= — q -^-VolyLog\?>, -z q ] 



(7) 
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FIG. 1: Behavior of the Energy density for pure gauge theory as a function of temperature 
Recall that the corresponding ideal value e 7 /T 4 reads: 
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The behavior of the energy density (e/T 4 ) as a function of temperature (T/T c ) is shown 
in Fig.l for pure gauge theory and for full QCD with iV^ = 2, 3 in Fig. 2. From these plots 
it is clear that e/T 4 approaches the ideal value only asymptotically. 



B. The entropy-density 



We compute the entropy density for the interacting pure QCD as well the interacting 
quark-gloun plasma. This is again a straight forward exercise since we have the equilibrium 
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FIG. 2: Behavior of the Energy density in 2- and 3-flavor QCD as a function of temperature 




FIG. 3: Behavior of the Energy density and entropy density as a function of temperature in Lattice 
QCD. 

distribution function for the quasi-partons already in hand. The entropy density in terms 
of the grand canonical partition function reads: 

S = lcMTln(Z g )] + icMTln(Z M )] 
ln(Z g ) = -Vv g J 01n(l-^exp(-^)) 
ln(Z M ) = ln(Z q ) 

d 3 p 



/dp 
— \n(l + z q exp(-/3p)). 



(9) 



The expression for the gluonic and quark contributions to the entropy density are obtained 



as, 
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FIG. 4: Behavior of the entropy density for pure gauge theory as a function of temperature 
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FIG. 5: Behavior of the entropy density in 2- and 3-navor QCD as a function of temperature 
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(10) 

The total entropy can be obtained by adding the gluon and quark contributions, S = 
S g + S q . We have plotted the dimensionless quantity S/T 3 for pure gauge theory and full 
QCD (Nf = 2,3), as a function of T/T c for EOS1 and EOS2. These are shown in Figs. (4) 
and (5) respectively. We shall utilize the expression for entropy density displayed in EgfTOl 
for determining the viscosity to entropy density ratio. The corresponding ideal value for 
S/T 3 is given by, 
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C. Comparison with lattice results 



We compare the thermodynamic observables determined by employing this model with 
the lattice results obtained by KarschQ]. The lattice results 42] are shown in Figs 3, and 



our results are displayed in Figs 1,2 and 4-5. For 2- and 3-flavor lattice results, we consider 
Fig. 14 of Ref. Q|. The agreement is overall good for EOS2, particularly beyond 2.5T C ; the 
agreement is merely qualitative for EOS1. These findings are consistent with our earlier 
result {20! that EOS2 predictions for the temperature dependence of the Debye mass and 
the dissociation temperatures of heavy quark ground states are in broad agreement with the 
lattice values. 



IV. VISCOSITY OF INTERACTING QUARK-GLUON PLASMA 

A. A brief review 

The determination of viscosity is not as straight forward an exercise as the determination 
of the thermodynamic observables. For, it requires modeling beyond the equilibrium prop- 
erties, in terms of the collision terms and other transport parameters, and also the nature 
of the perturbation to the equilibrium distribution. We note here that for sQGP under 
consideration, the collision term is by no means easily determined since perturbative results 
involving lowest order contributions are by no means guaranteed to be reliable. A reliable 
non-perturbative collision term is even harder to obtain. We need to adopt methods that 
go beyond the determination of properties such as the energy density or the specific heat. 



theor 



There are two ways to compute the transport parameters for QGP:(i) from quantum field 
using the Kubo formula^, [j^] or (ii) from the semi-classical transport theory {13I. 



14|,l32|, 



331 ] . To model the QGP produced in heavy ion collisions employing semi-classical 
transport theory, one needs to employ the Vlasov term which incorporates the dynamics 
of non-abelian color charges. We also need a reliable collision term, which is, as we have 
pointed out, difficult to determine. A collision term which has recently been computed 



by Arnold, Moore and Yaffe 



32] by considering binary collisions at the tree level in the 



lowest order. The collision term so obtained is then utilized to estimate the shear viscosity 
for QGpflE]. 

Further, Asakawa et. al [l_4j have included the Vlasov term for an ensemble 
of turbulent color fields, and determined the anomalous shear viscosity;this determination 



does not require a collision term. This work generalizes the earlier work of Dupree 34| by 
including the nonabelian dynamics. Both these works are based on methods describe in 
33]. As mentioned earlier, studies based on AdS/CFT correspondence j^, 35] employ the 
Kubo formula and predict a lower bound of (l/An) for 7]/S for plasmas whose dynamics are 
governed by class of strongly coupled gaug e theories. This speculation is supported in some 



13 



16]. 



studies based on transport theory 

In a recent work [La. it has been argued that one does not need to treat QGP as a 
strongly coupled plasma to understand low viscosity. The authors argue that the anomalous 
transport processes in the rapidly expanding QGP are actually responsible for the very low 
value of the shear viscosity, and not the binary collisions. On the other hand in a very recent 
work, Xu and Greiner[l^] argue that the the reason for a small value of rj/S is mainly due 
to gluon bremsstrahlung contributions to the collision term. 

In the present paper, we adopt the approach of Asakawa et a/fl^] and determine the con- 
tribution of both the collisional and the anomalous parts to the ratio f]/S. As in the earlier 
study 15] , we find that the anomalous part dominates over the collisional contribution. 



B. Determination of r/ 

In this Section, we determine the viscosity of a rapidly expanding interacting QGP, as 
described by EOS1 and EOS2, and as represented by equivalent equilibrium distribution 
unctions (section 2). Our procedure involves replacing the ideal gas distributions used in 



141 1. by the ones obtained by us for EOS1 and EOS2. We find that all the assumptions 
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made in Refs. [141. 1151 ] will be applicable in the present case. 



Let us firs t briefly outline the standard procedure of determining viscosity in transport 
33] . The shear viscosity of QGP in terms of parton occupation numbers can be 



theory [1J. 

obtained by comparing the microscopic definition of the stress tensor with the macroscopic 
definition of the viscous stress tensor. The microscopic definition of the stress tensor in 
terms of the distribution function is as follows: 

Tlk = J J^E p PiPkm ^ (12) 
On the other hand the macroscopic expression for the viscous stress tensor reads: 

T ik = P5 ik + euiu k - 2r](Vu) ik - (5 ik V ■ u, (13) 



where rj is the shear viscosity, ( is the bulk viscosity and ( Vm)^ is the traceless, symmetrized 
velocity gradient, 



(V«)y = -(ViUj + VjUi) - -SijV ■ u. (14) 
To obtain the shear viscosity, we write the distribution function as 

r) = z g/f eM-PE p + h{p^)Tl (15) 
Assuming that fi(p,r) is a small perturbation to the equilibrium distribution, we expand 

f(p,r) and keep the linear order term in f\; the following form of the distribution function 

is thus obtained: 



f(p,r) = fo{p) + &f(p,r) 

= /o(p)[l + A(p,r)(l±/o(p)], (16) 

where f = jlq { see EqfJ]). 

The important point to be noted is that while deriving the transport coefficients, one as- 
sumes a slow variation of the particle distribution so that the deviation from the equilibrium 
distribution is homogeneous in space and proportional to the gradients of the equilibrium 
parameters. Employing the standard approach Q, Q, we write 

m^ = -^ PtPj (vu) tv (17) 

where the dimensionless function A(p) measures the deviation from the equilibrium config- 
uration. Since t] is a Lorentz scalar, it may be evaluated conveniently in the local rest frame. 
For a boost invariant longitudinal flow, (V«)y = -^diag(—l, —1, 2) in the local rest frame, 
and and fi(p) takes the form 



A(p) / 2 V 



A® = -3#Hrf-T)- < 18 > 



where r is the proper time(r = \Jt 2 — z 2 ). 
The expression for t] is then obtained as 



V = ^- f ^^A(p)^ (19) 

V 15 J WE* KP) dE p ' [ J) 



entirely in terms of the unknown function A(p). E p is the particle energy. 
We adopt the ansatz in jbfl and take the form of A(p) to be 

A(p)=A\p\/T-A = A g ,A q . 



(20) 



To appreciate the ansatz better, we note that the dispersion relation in EqJl]) gets modi- 
fied, in the presence of the perturbation fi as given by 



E eff(p) = E p 

The velocity is given by (v p = dpE e ff(p)) 

v p = p 



pA 



E p T 3 t 
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From this expression for v p , it is clear that the perturbation leads to different velocities in 
transverse and longitudinal directions. This introduces manifest anisotropy in the system. 
We determine A(p) by the variational method by minimizing the linearized tr ansp ort 



n 



equation |33| with a Vlasov term and a collision term computed by Arnold et. aZ[32j. 



The factor A is yet undetermined. To fix its value, we minimize the quadratic functional 
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(23) 



where the first term gives the drift, the second represents the diffusive Vlasov dynamics, and 
the last term is the collision integral. The expression in parenthesis is just the transport 
equation satisfied by h(p) after averaging over the color fields JjJ|. The equilibrium dis- 
tribution functions modify the results obtained in [ljj] by rendering the coefficients A q , A g 
dependent on temperature and the coupling constant. After performing the momentum in- 
tegrals, the drift, diffusive and the collisional terms in the quadratic functional Eqj23l acquire 
finally the form 
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In = -PolyLog[n, -z' 1 ]. 



The function PolyLog[n,a\ has the series representation 

OO I, 



(25) 



PolyLog[n,a} = Y J -^ l - W 

k=0 



These expressions reduce to the ones obtained in [l^], if we put z q , z g = 1, corresponding 
to ideal quark and gluon distributions. 

C. The anomalous and collisional viscosities 

Let us now turn our attention to determine the analytic expressions for anomalous and 
collisional contributions to the shear viscosity, which are determined respectively by the 
diffusive Vlasov and the collision terms in EqfSH ■ To determine either of them, one utilizes 
Eqll9] along with EqEUby following exactly the path taken in Reffl^. 

By inserting Eq{19] in Eqj20] and performing the momentum integration, we obtain the 
following expression for viscosity r\: 

v = v^~" [ {N * ~ 1)71 Ag + NcNfI! Aq \ ■ (27) 

The minimization of the functional I^[/i](Eq j2"4l leads to the following matrix equation 
for the column vector A = (A g , A q ): 



(a A + ac)A = f, (28) 



where, the column vector r and the matrices and ac are given by 

_ 32 / {Nl - 1)1! 
r = I 

3n2 \ N c N f II 
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where C g = (TV 2 — l)g 4 log(g -1 ). This leads to the following expressions for the anomalous 
and collisional contributions to the shear viscosity, 

3 



3 



Vc = ^p- 3 r-(ac 1 )-r (32) 
By employing the additivity of rates, the expression for the total viscosity is obtained as 

Q 

± = ± + J-. (33, 

It is clear from Eqs 1291(301 I3T1 and [32] that r) c ~ g^wI/gT- I n the wea k coupling limit 
(g << 1), at which the hot EOS are also valid, the total viscosity in Eqj33]will be dominated 
by the anomalous component. Therefore for a weakly coupled QGP rj ~ tja- Hence we 
confine our attention to anomalous shear viscosity and study it's behavior with temperature. 
The individual expressions for the gluon and quark contribution to the anomalous viscosity 
from EqJHH are obtained to be 
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(34) 



The total anomalous shear viscosity is obtained by summing up these two contributions, 
(i.e, = r] 9 A -\- r) A ). We note that the above expressions are valid for a purely magnetic 
plasma. For the case when both chromo-electric and chromo -magnetic fields are present in 
the turbulent phase, and all their components are of equal size, the expressions for viscosity 
can be obtained simply by the replacement < B 2 > r m — > |(< S 2 > + < B 2 >)r m fl^|. 
Accordingly, we rewrite EqJ3H as 

lo/r 6 (n 2 - i) 2 (/f) 2 
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We pause to compare the viscosities with their ideal values. Recall that the contribution 
from ideal distribution functions are obtained by setting z g = 1 and z q = 1 in EqsJ34landl35l 
as expected they match with the expressions of Asakawa et. a/|l^|. Thus, the expressions 
for the relative viscosities (anomalous) are read off as 
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The behavior of lZ q and 1Z q as a function of temperature for EOS1 and EOS2 are shown in 




FIG. 6: Behavior of the relative viscosity for pure gauge theory as a function of temperature 



Figs. (7) and (8). Clearly incorporation of interaction effects in the EOS further reduces the 
viscosities. While the gluon viscosity can reduce upto ~ 30%, the fall in the quark viscosity 
can be as steep as 80%, indicating a more ideal fluid like behaviour. 
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FIG. 7: Behavior of the quark contribution to the viscosity in 2- and 3- flavor QCD as a function 
of temperature 

D. Viscosity to entropy density ratio 



A determination of the absolute values of viscosity requires further a knowledge of the 
quenching parameter <7,fj[l(j|, which is defined as the rate of growth of the transverse mo- 
mentum fluctuation of a fast parton in an ensemble of turbulent color fields In turn, 
q R is given by 



<1r 



■{E 2 + B 2 )r n 



(37) 



3(iV* - 1) 

in terms of the total energy density and an appropriate relaxation time r [14] . 

One can combine this with the expression for the anomalous viscosity of gluons and obtain 
the relation as follows, 

v 9 A(z 9 ) = ^W-i)& 6i 4r- (38) 



Estimates for qn are available 



or the gluonic case only. In this case, qn is estimated from 



the data by various approaches 36]. Studies within the framework of the twist expansion by 
: itting the experimental data on hadron suppression in the most central Au-Au collisions 



23, 



yield values in the range 1 — 2GeV 2 / fm for the gluon quenching parameter. On 
the other hand, an eikonal approach (39), l^j estimates it to be roughly ten times larger than 
the twist estimates, in the range 10 — 30GeV 2 / fm. The expression for the gluonic rj/S is 



obtained as 



(/f) 2 



V = 20T e 3 (A r e 2 — 1) T 3 

S ~ 3 ^ T c I*[v g (4II-\n(z g )l!) + '' 



(39) 

'3ryv"4 — \~yj-6J 1 T 4 J 

It is clear that the estimates for rj/S inherit the uncertainty in g^, upto an order of mag- 
nitude. For purposes of concreteness, we choose the QCD transition temperature (T c ) to be 



270Afev[2g]. We have plotted EqJH]as function of temperature for E0S1 and E0S2 in Fig.8 
and 9, with respective values £r = 1, lOGeV 2 / fm. We see that in the latter case, the ratio 
can fall significantly below the AdS/CFT bound ^- ~ 0.08 even at 3T C , but it may not be 
reliable since the large value of (qn = lOGeV 2 / fm) which we have employed may not be 
accomodated within weak coupling framework 16] which we consider in the present paper. 
On the other hand, in the former case = lGeV 2 / fm the value of rj/S does not violate 
the AdS/CFT bound, although it is quite close to it near 2T C . It appears that the violation 
of bound, which can occur at q^ > \lGeV 2 / fm will be marginal near 2T C . 

As expected the ratio increases with increasing temperature. Interestingly, unlike other 
thermodynamic variables, the ratio is not sensitive to the EOS employed. 
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FIG. 8: Viscosity to entropy density ratio for pure gauge theory for or = IGeV 2 / fm 

We now establish the connection between our results and that of Asakawa et al [jjj] . The 
expression for the ratio rj/S in Eqj39] reduces to that of Asakawa et al if we set z g = 1 and 
employ the ansatz, A(p) = A g / q p/T for the anisotropy parameter 43|. The corresponding 
expression in this limit reads, 



V 20T!(N 2 

s 



1),2\b C(5) ! 



(40) 



u g q R y T c > 4C(4) 2 ' 

Let us normalize the viscosity to entropy ratios for EOS1 and EOS2 wrt the ideal values. 
We have plotted the relative ratios, which we denote by 1Z V , as a function of temperature, 
in Fig. 10, which shows the effects of interactions in rj/S. From Fig.10, we see that (7^) is 
less than unity, approaching the ideal value asymptotically. Interestingly, the EOS2 values 
are closer to the ideal case, differing by about 3% near 2T C . 

Finally, note that the expression for 7]/S in EqfJU] is identical to the expression used in 
except for a numerical factor of 0(1). This discrepancy arises because we consider both 



;he diffusive Vlasov and the collision terms in the transport equation, while the analysis of 
Ji| neglects the collision term. 





T/T c 



FIG. 10: Behavior of TZ V as a function of temperature in the case of pure gauge theory for EOS1 
and EOS2. Note that 1Z V scales with T/T c . 



V. CONCLUSIONS AND OUTLOOK 



In conclusion, we find that hot QCD EOS corresponding to interactions of 0(g 5 ) and 
0(g 6 ln(l/g) can significantly impact the values of the thermodynamic observables such as 
the energy density. The viscosity and the ratio r]/S, which we have studied as functions 
of temperature, get reduced by approximately 7% for EOS1 and 4% for EOS2 near 2T C in 
contrast to their ideal counterparts. We found that the value of rj/S for q^ = lGeV 2 / fm 
near 2T C is closer to the lower bound 1/47T placed on t]/S by AdS/CFT studies. Further, 
the choice qn ~ 10GeV 2 /fm is difficult to accommodate within the weak perturbative 



framework and hence the violation of the AdS/CFT bound may not represent the factual 
situation. The choice ~ 2GeV 2 / fm does lead to a violation near 2T C , but only marginally 
so. In short, the findings in the present work strengthen the near perfect fluid picture of 
the hot and dense matter created in relativistic heavy ion collisions. This analysis has 
been rendered possible because of the mapping of interacting partons to non-interacting 



quasi partons with effective fugacities 19|, |20j. While the present study points definitively 



to the importance of interaction effects, it is by no means complete, because of inherent 
uncertainties in the estimates of the gluonic quenching parameter, and an absence of the 
knowledge of the quenching parameter for the quarks. The EOSs which we study are also 
perturbative. It should be of great interest to employ the lattice EOS^, \^\. This will be 
taken up separately. 

Acknowledgment : 

We are thankful to Frithjof Karsch for providing us lattice data for pure gauge theory 
and Anton Rebhan for useful comments and suggestions on fixing 0(g 6 ) contributions at 
0(g 6 \n(l/g)). VC thanks Raman Research Institute for hospitality where the manuscript 
was finalized and CSIR, New Delhi(India) for financial support. 



[1] STAR Collaboration , Nucl. Phys. A 757, 102 (2005), PHENIX Collaboration , Nucl. Phys. 
A 757, 184 (2005), PHOBOS Collaboration, Nucl. Phys. A 757, 28 (2005), BRAHMS Col- 
laboration , Nucl. Phys. A 757, 01 (2005). 



[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 
[10] 



PKovtun, D.T.Son, A.O.Starinets, Phys. Rev. Lett.94 , 111601 (2005). 
R. Baier and P. Romatschke, |arXiv:nucl-th/0610108| 
Drescher et al, Phys. Rev. C 76, 024905 (2007). 



Frithjof Karsch, Lect. Notes Phys. 583, 209 (2002) ( a rXiv : hep-lat/0106019[ ) . 

E. V. Shuryak, Nucl.Phys. A 774 387 (2006), [hep"-ph/0608 177] Boris A. Gelman, Edward V. 

Shuryak and Ismail Zahed, Phys. Rev. C 74, 044908 (2007) (|arXiv: 0704 .35531 (nucl-th)). 

H. B. Meyer, Phys. Rev. D 76, 10171 (2007) (jarXiv : 0704 . 180H Chep-lat ) . 

A. Nakamura and Sunao Sakai, Phys. Rev. Lett. 94, 072305 (2005). 

Adare et al, Phys. Rev. Lett. 98, 172301 (2007). 

Sean Gavin and Mohamed Abdel-Aziz, Phys. Rev. Lett. 97, 162302 (2006). 



[11] STAR Collaboration , Nucl. Phys. A 757, 102 (2005). 

[12] Lacey et al, Phys. Rev. Lett. 98, 092301 (2007). 

[13] Zhe Xu and Carsten Greiner, Phys. Rev. Lett. 100, 172301 (2008). 

[14] Masayuki Asakawa, Steffen A. Bass and Berndt Muller, Prog. Theor. Phys. 116, 725 



(2007) QarXiv : hep-ph/0608270) . 



[15] Masayuki Asakawa, Steffen A. Bass and Berndt Muller, Phys. Rev. Lett. 96, 252301 (2006). 
[16] Abhijit Majumdar, Berndt Muller and Xin-Nian Wang, 



Phys. Rev. Lett. 99, 192301 (2007) ( arXiv : hep-ph/0703 082 ) . 



[17] P. Arnold and Chengxing Zhai, Phys. Rev. D 50, 7603 (1994); Phys. Rev. D 51, 1906(1995); 

Chengxing Zhai and B. Kastening, Phys. Rev. D 52, 7232(1995). 
[18] K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Phys. Rev. D 67, 105008 (2003). 
[19] Vinod Chandra, Ravindra Kumar and V. Ravishankar, Phys. Rev. C 76, 054909 

(2007 )arXiv : 0705 . 2690i rmucl-thl . 
[20] Vinod Chandra , Akhilesh Ranjan and V. Ravishankar, arXiv:0801 . 1286 [hep-ph] . 



[21] Anton Rebhan hep-ph/0504023 



[22] S. Haung and M. Lissia, Nucl. Phys. B 438, 54 (1995). 

[23] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983). 
[24] E. Braaten and A. Neito, Phys. Rev. D 53, 3421 (1996)(|nep z ph/9510408"l). 



[25] A. Vuorinen, arXiv:hep-ph/0402242. 



[26] J.-P. Blaizot, E. Iancu, A. Rebhan, Phys. Rev. D 68, 025011 (2003). 
[27] S. Gupta, Phys. Rev. D 64, 034507 (2001)(hep-lat/oolooll). 



[28] Olaf Kaczmarek and Felix Zantow, Pos(LaT 2005)177( [arXiv : hep-lat/0510 093 

[29] Akhilesh Ranjan and V. Ravishankar, larXiv : 0708 . 3697l (nucl-th) . 

[30] Helmut Satz, Nucl. Phys. A 783, 249 (2007). 

[31] W. M. Alberico, A. Beraudo, A. De Pace and A. Molinari, Phys. Rev. D 75, 074009 (2007). 

[32] P. Arnold, G. D. Moore and L.G. Yaffe, JHEP 0011, 001 (2000). 

[33] E M Lifshitz and L P Pitaevskii, Physical Kinetics (Landau and Lifshitz; Volume 10) Pergamon 
Press. 

[34] T. H. Dupree, Phys. Fluids 9, 1773 (1966), Phys. Fluids 11, 2680 (1968). 

[35] Buchel A, Liu J T and Strainets A O, Nucl. Phys. B 707, 56 (2005). 

[36] A. Majumder, arXiv:nuc l-th/0702066, 



[37] H. Zhang, J. F. Owens, E. Wang and X. N. Wang, |arXiv : nucl-th/07010451 
[38] A. Majumder, C. Nonaka and S. A. Bass, larXiv :nucl-th/0703019l 

[39] R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schi, Nucl. Phys. B 483, 291 
(1997). 

[40] N. Armesto, L. Cunqueiro, C. A. Salgado, W.-C. Xiang. larXiv : 0710 . 30T3l (hep-ph) 

[41] Frithjof Karsch, larXiv : 0711 . 0661li rhep-latl ; M. Cheng et al, Phys. Rev. D 77, 014511 

(2008) (jarXiv : 07l0 . 0354i rhep-latl ) . 
[42] We thank Frithjof Karsch for providing us the pure lattice gauge theory data. 
[43] In writing EqHUJ we have employed the following results: VolyLog[4, z g = 1] = ("(4) = 7r 4 /90 
and TolyLog[5, Zg = 1] = £(5) and that S 1 = 4^C(4)T 3 /vr 2 . Eqg0]follows from from Eqs.(6.32 



and 6.33) in Ref. 



14] by employig the form of given above. 



